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Chandra observations. 
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^ [ We present new interpretation of Ton S180 spectrum obtained by Chandra 

I Spectrometer (Low Energy Transmission Grating). Several narrow absorption 

I lines and a few emission disk lines have been successfully fitted to the data. We 

^ ' have not found any significant edges accompanying line emission. We propose 

' the interpretation of narrow lines consistent with the paper recently written by 

Krolik (2002), where warm absorber is strongly inhomogeneous. Such situation 
is possible in so called multi-phase medium, where regions with different ion- 
ization states, densities and temperatures may coexist in thermal equilibrium 
under constant pressure. We illustrate this scenario with theoretical spectra of 
radiation transfered through a stratified cloud with constant pressure (instead 
of constant density) computed by code titan in plane parallel approximation. 
Detected spectral features are faint and their presence do not alter the broad 
band continuum. We model the broad band continuum of Ton S180 assuming 
an irradiated accretion disk with a dissipative warm skin. The set of parameters 
appropriate for the data cannot be determined uniquely but models with low 
values of the black hole mass have too hot and radially extended warm skin to 
explain the formation of soft X-ray disk lines seen in the data. 



- 2 - 



Subject headings: X-ray: galaxies - galaxies: Seyfert -galaxies: individual: Ton 
S180 



1. Introduction 

Narrow Line Seyfert 1 galaxies (Osterbrock & Pogge 1985) form a special class of active 
galactic nuclei. Narrowness of the optical lines frequently correlates with the steepness of 
the soft X-ray spectra (Boiler, Brandt & Fink 1996) hinting for a connection of this class 
with soft state Galactic sources. The overall spectra are sometimes well described by a 
disk emission superimposed on a broad IR-X-ray power law component (Puchnarewicz et al. 

2001) but in most cases the soft X-ray slope seems to be too shallow for an extension of a disk 
black body emission. Such a steep soft X-ray spectrum of a roughly power law shape may 
form either as a result of Comptonization (e.g. Czerny & Elvis 1987) by a moderately hot 
material (in comparison with the very hot plasma responsible for the hard X-ray component 
in normal Seyfert 1 galaxies) or as a result of reflection by a strongly ionized material (e.g. 
Czerny & Zycki 1994, Czerny & Dumont 1998). The established presence of warm absorber 
seen in many of those objects (CoUinge et al. 2001, Kastra et al. 2002, Kaspi et al. 2002, 
Kinkhabwala et al. 2002, Leighly et al. 2002, Yaqoob et al. 2002) complicates the analysis. 
However, with the extremely good spectral resolution of Chandra and XMM the issue may 
be addressed through the careful analysis of the existing spectral features. 

Ton S180 {z — 0.06198, corresponding to 286 Mpc for Hubble constant 65 km s^^ 
Mpc~^) is an extreme case of a NLSl with its low value of the FWHM of Hf^ line (~ 900 
km s~^; Wisotzki et al. 1995) and steep soft X-ray spectrum of photon index F = 2.68 
(Comastri et al. 1998) The Galactic absorption toward the object is low {N^ = 1.52 x 10^° 
cm~^; Stark et al. 1992). The source was extensively observed in all energy bands. The 
broad band spectral energy distribution (SED) determined by Turner et al. (2002) indicates 
that most of the energy is emitted in the band 10 — 100 eV. Furthermore, HST data show 
typical emission lines as expected for a Seyfert 1 galaxy, and the lack of absorption features. 
High resolution spectra obtained by FUSE reveal UV absorption by Ovi and the absence of 
neutral hydrogen absorption, indicating high ionization state for the absorbing gas. 

High resolution X-ray spectra of this object were recently obtained with XMM-Newton 
and Chandra sateUites. XMM-Newton spectrum was steep and featureless (Vaughan et al. 

2002) . Chandra data also appeared predominantly featureless and well represented by a 
single power law with some soft excess (Turner et al. 2001, Turner et al. 2002). The shape 
of the soft excess was fitted by a DISKBB model, however plotted residua hinted for the 
possible presence of sharp spectral features (Fig. 1 Turner et al. 2001). 
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Therefore, in the present paper we reanalyze Chandra data in some detail (Sec. 2 and 
3), using a new response matrix acisleglD1999-07-22rmfN0002.fits and including any effects 
related to the contamination layer on ACIS-S. We note that the observation of Ton S180 
was performed on day 143 of the Chandra mission when the effects of the contaminant were 
much smaller than currently observed. We detected several weak narrow absorption lines 
due to the material located at the line of sight to the source. Also, several individual disk 
lines are fitted with equivalent widths in rough agreement with the predictions of the model 
of illuminated accretion disk spectrum presented by Rozanska et al. (2002 a,b). 

In Sec. 3.2.1 we discuss existing models (Kaastra ct al. 2002, Kinkhabwala ct al. 2002 
and others) for fitting absorption by warm, partially ionized gas present in many AGN. We 
point out, after Krohk (2001), that the warm absorber may not be the slab of constant density 
matter, but it is rather the distribution of regions with different densities and temperatures. 
Such a multi-phase mixture of gas can exist under constant pressure, as a result of thermal 
instabilities (Krohk, McKee & Tarter 1981). 

Till now, the multi-phase models of warm absorber were not constructed, because pop- 
ular photoionization codes (xSTAR, cloudy), compute radiative transfer through the con- 
stant density matter. In Sec. 4, we present new set of models computed by the code titan 
(Dumont, Abrassart & Collin 2000), which solves radiative transfer through the matter in 
non-LTE approach. Introducing assumption of constant pressure (not density) of warm ab- 
sorber, we obtain spectra with many absorption lines from the gas with strong temperature 
and density gradient. Our code is suitable for any optical thickness of the absorber, never- 
theless, we consider regions with total column density of order of lO^^cm^^, as suggested by 
Krolik (2001) to be a typical property of warm absorber. 

In Sec. 5 we construct two global disk corona models trying to predict the slope of ob- 
served continuum. The model with higher black hole mass (M = IO^Mq) is more promising 
in the sense that emission lines from the disk are expected in the spectra, since hot skin does 
not cover disk completely. 

Sec. 6 is devoted to the discussion, while conclusions are listed in Sec. 7. 

2. Chandra observations and data reduction 

Ton S180 was observed with the Chandra satellite in December 14, 1999, for about 80 
ksec in Low- Energy Transmission Grating in ACIS-S/LETG configuration, as described by 
Turner et al. (2001) and Turner et al. (2002). The resolution of LETC corresponds to 
an FWHM pa 1200i?fcey km s~^ for the range of energies from 0.1 up to 6 keV. The source 
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luminosity at this epoch (2-10 keV unabsorbed flux equal 4.6 x 10~^^ erg s"^ cm~^ and 
0.3-10 keV unabsorbed flux equal 2.2 x 10~^^ erg cm~^ ) was comparable to that during 
the iJe^^o-^/lX observation (2-10 keV: 4.2 x 10"^^ erg s~^ cm~^; Comastri et al. 1998) and 
XMM-Newton observation (0.3 - 10 kcV: 2.2 x 10"^^ erg s~^ cm'^; Vaughan et al. 2002) and 
somewhat fainter than during the 12 days ASCA monitoring (2-10 keV: 6.5 x 10~^^ erg s~^ 
cm~^; Romano et al. 2002). 

We obtained the data from the Chandra archive. The data were processed with the stan- 
dard CXC data processing pipeline ASCDSVER 6.7.0. We analyzed the data using ClAO 
Version 2.2. and follow Grating Analysis Thread (CXC web page) to extract the spectrum 
and create appropriate effective area files (ARF files). We combined both sides of the first or- 
der spectrum. We used calibration data base version CALDB v.2.14, in particular we apphed 
acisD 1997-04- 17qeN0003.fits, quantum efficiency file and acisleglD1999-07-22rmfN0002.fits 
as response matrix. This formally allows us to study the entire energy range between 0.12 
till 12.2 keV, although calibration uncertainties are higher below 0.3 keV, while the effective 
area drops significantly above ~ 4 keV. After removing bad pixels and rebinning the data to 
achieve at least 20 counts per channel we examined the energy between 0.18 and 6 keV. All 
spectral fits have been performed with the XSPEC 9.0 fitting package. 

Since the bcgining of the Chandra mission ACIS-S detector has been accumulating a 
layer of a contaminant which affects any measurements at low energies (below ~ 0.8 keV). 
Detailed effects of the contaminant arc still being investigated by the Chandra X-ray Center 
^, however, it is clear that the overall degradation of ACIS quantum efficiency and presence 
of additional absorption features is time dependent. Ton S180 was observed at the day 
143 of the Chandra mission and based on the currently available calibration data we expect 
that the instrument shows maximum of ~40-15 % degradation within 0.28-0.6 keV range 
(i.e. maximum ~ 40 % in 0.28 keV and ~ 15 % in 0.6 keV). In the following we identify 
the possible absorption features due to the absorber in Ton SI 80 using the entire available 
energy range. We note here that some of the features at energies close to the edges of the 
contaminant (C,N,0,F K-edge) may be affected. We use the available model, acisabs, to 
include the contamination effects in the analysis. 



^http: / / asc.harvard.edu/cal/Acis/CaLprods/qeDeg/index.html 
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3. Spectral analysis 
3.1. Continuum 

In the first step, we fit the Nh value of Galactic absorption together with a broken power 
law continuum. Our best fit value of NH{Gai.) — 1-83 ± 0.15 x 10^° cm~^ is slightly larger 
than Galactic absorption NH{Gai.) — 1-5 ± 0.5 x 10^° cm~^ determined by Dickey & Lockman 
(1990) and used by Comastri et al. (1998) or than NH{Gai.) = 1-52 x lO^*' cm~^ determined 
by Stark ct al. (1992) and used by Turner et al. (1998, 2001). All fitted parameters for 
"plain continuum" are listed in Table 1. 

In the next step we fit broken power law with value of Galactic absorption fixed at 
NH{Gai.) — 1-52 X 10^'^ cm^^ (see Table 1 second line). It results in slightly higher y^, but the 
power law index below 2 keV is Fi = 2.69 ±0.02, and agrees with the best fit slope V = 2.68 
presented by Comastri et al. 1998 for Beppo SAX data. Lower value obtained by Turner ct 
al. (2001 and 2002), equal to 2.44 ± 0.04 resulted from the use of additional component to 
model the soft X-ray excess. 

The energy break obtained by us in both models (about 1.7 ± 0.14 keV) is somewhat 
lower than the value determined by Comastri et al. 1998 {Ebreak = 2.5), but the spectral 
slope above this energy break, r2 = 2.26 ± 0.07 agrees with results Comastri et al. 1998. 

Since both broken power-law indices are very similar to those obtained from other 
satellites we adopt second model (fixed Galactic absorption) for further analysis, leaving 
= 1197 per 1682 d.o.f. (reduced = 0.7105). 

After fitting continuum, there is still some structure in the residuals of the soft X-ray 
spectrum, suggesting the presence of absorbing matter along the line of sight toward the 
source. Therefore, we have tried to determine column density of this gas in the simplest way, 
by fitting Ovil, Ovill, and Cvi edges often seen in Seyfert galaxies (Mrk 509 Yaqoob et al. 
2002, NGC 3783 Kaspi et al. 2002). Surprisingly, such edges appeared to be very weak in 
the data of Ton S180. We were able only to give upper limits for rovii ^oviii, and rcvi of order 
of 10~^, which practically indicates the lack of those edges. 

3.2. Spectral features 

3.2.1. Application of standard warm absorber 

X-ray absorption lines appear to be common feature in many AGN. The widely accept- 
able explanation is the existence of warm, ionized gas which is illuminated by the energetic 
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raxliation from nucleus. Trying to find properties of warm absorber we need proper tools to 
fit those data. Till now, models of absorption by the partially ionized plasma available in 
XSPEC (absori, photo,) do assume slab with constant density. 

In ABSORI, the radiative transfer of lines is not computed at all, only photoionization 
cross sections are determined and absorption "on the spot" is computed Fout — Finexp^—Tune) 
(Doneetal. 1992). Constructed by Kinkhabawa et al. (2002) local photo model (http: //xmm.- 
astro.columbia.edu/research.html) neglects photo-electric absorption, but allows for individ- 
ual line absorption, treating each ionic column density separately with an initially unab- 
sorbcd power law. This method produces consistent results only when absorbing matter is 
optically thin, and transmitted spectra do not display photo-electric edges. Furthermore, 
they have assumed that plasma is optically thin to reemitted photons. 

In the aim to find evidences for warm absorber in Ton S180, we have applied ABSORI to 
our data. Allowing for three parameters to be fitted in this model, decreases by 38 (now 
reduced ^ = 0.747). The fit indicates that the medium has column density 7.34lo;22 ^ 10^^ 
cm~^ and ionization parameter (^ = 4178l^32|) for the fixed default temperature Tabs — 
3 X 10^ K. This model, however, leaves considerable residual trends and does not allow for 
a deeper insight into the ionization state of the absorbing material. It is not surprising 
since the model is based on simple computations of photoionization cross sections without 
accurate transfer through the medium. 

We have apphed photo to our analysis. The model again indicates very low column 
densities of each ion, about 10^^ cm~^ for Cvi, Ovil, and Ovill. We conclude that most 
probably, the warm absorber in Ton S180 is more complicated than optically thin constant 
density slab. It may be the stratified density gas, where accurate radiative transfer, including 
all photoionization processes, should be calculate to model transmitted spectrum. Such 
models are not available for fitting the data yet, so two or three separate absorbers with 
different constant densities and ionization states are used together to explain spectral features 
in some AGN (Kaastra et al. 2002 , Kaspi et al. 2002). This motivates us to present in 
Sec. 4 the model of the multi-phase gas with stratified density and temperature, where the 
radiative transfer of incident power law is carefully computed. The model is not ready for 
fitting yet. but gives interesting results based on equivalent widths of various lines. Observed 
line properties in the Ton S180 data are determined in the section below. 

More advanced modeling was done with the code XSTAR (developed by Kallman & 
Krolik 1999) which allows to estimate the total column densities of absorbed gas. Models 

are computed assuming constant density slab. Such models are very useful for reconstruction 
of the fiat part of curve of growth (Yaqoob et al. 2002), and for subsequent estimation of 
turbulent velocities of the absorber. Nevertheless, at the last stage of analysis separate 
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Gaussian components are fitted to the data. 

3.2.2. Absorption lines 

We have successfully fitted several absorption lines using standard Gaussian model for 
each line separately (see Figs. 1, 2, 3, 4, 5, 6, and 7. Gaussian fit requires three free 
pcirameters, energy of the line centroid, Gaussian width and normalizations, 
which has to be negative in case of absorption Unes. We include in the fit only those 
features which are significant. All absorption lines are presented in Tabic 2. For each line 
we list its central energy in keVs, equivalent width (EW) in eVs, the decrease of x^, best 
identification with observed and rest frame energy, and with possible blends. 

Some of those fines, especially the softest ones, may be caused by dust presented on 
Chandra mirrors. This should show up a common feature in all objects even those with 
different rcdshifts. To check whether lines are real or instrumental, we compare observations 
of Ton S180 with Chandra data of quasar 3c273 observed two weeks latter with the same 
configuration as Ton S180. Exemplary residuals (for soft range between 0.2-0.24 kcV) for 
both objects are shown in Fig. 8. In both cases local power-law was fitted with fixed value of 
Galactic absorption appropriate for each object. The same features around energy 0.22 keV 
are present as well in 3c273 with redshift 0.156 as in Ton S180. We repeated such analysis 
in the whole energy range. All possible instrumental fines are labeled in Table 2. 

Some of the detected features are seen only in Ton SI 80 spectrum and therefore seem to 
be real, intrinsic to the source. For example, in Fig 9 we show that line Cvi with observed 
energy 0.4428 keV which is detected in Ton S180 data is not seen in 3c273. 

Lines are very faint. They usually have equivalent width less or about 2 eVs, and 
several of them are just marginal detections, as clearly seen from the line EW error, given 
at 90% confidence level (each error is computed for the one parameter of interest). 

Only two of them are stronger, but from theoretical predictions by analysing our modeled 
spectra from Sec 4, they are not blends. All possible blends according of our modeled specta 
obtained with spectral resolution 100 are listed in Table 2. There is therefore no major 
discrepancy between the Chandra results and the XMM observations (Vaughan et al. 2002) 
where only upper limits of order of 1 eV at energies 0.35 - 0.7 keV were provided. Most lines 
originate from hydrogen or helium like ions indicating that matter is highly ionized. There 
is an absence of lines from mildly ionized spices. 

Interestingly, we found two different transitions from the same element at the given 
ionization level, (Cvi(i_2) and Cvi(i_5) where i - j denotes number of levels from and to 
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which the transition proceeds). This is helpful in the aim to estimate the optical thickness 
of Cvi ions and then total optical thickness of the warm absorber. 

If the medium is optically thin in those lines (i.e. weak line case), we are on the linear 
part of the curve of growth, where equivalent width is proportional to the column density 
Ni of the absorbed ion: 

EW, = NiUj, (1) 

rrieC 

where fij is the oscillator strength transition and stimulated emission has been set equal 
to unity. Thus, for the given population of ions at the same ionization level, the ratio of 
equivalent widths is given by: 

(2) 



EWik fik 

Since the oscillator strengths are strictly atomic features independent from the radiation 
field and are calculated theoretically we use them to find the ratio 

On the other hand, the same ratio determined for Ton S180 (see Table 2) equals 0.62l|^^. 
Error is computed assuming most pessimistic case, when numerator has its maximal value 
EW + ^EW ^ and denominator has its minimal value EW — AEW. Although the error is 
huge, the value of the ratio determined from observation is still far from theoretical predic- 
tions for optically thin medium (Eq. 3). Our theoretical spectra (see Sec. 4) suggest that 
C(i_2) line is possibly a blend with Mgxil, Sixiv and most strongly with Cv. We have 
computed that real EW of the C(i_2) line can be even up to four times lower than the EW 
of the whole blend. But this fact even lowers the observed ratio of EWs of two carbon lines. 

The straight conclusion is that the emitting region is optically thick for those lines and 
the transfer in lines should be taken into account for proper modeling of the warm absorber. 
It is also possible that the identification of the lines is not correct or the line is a blend of 
the carbon line with another transition. Nevertheless, we further explore the possibility that 
the carbon line ratio is determined correctly. 

Therefore, the Equation 1 gives only a lower limit for the ion column density, equal 
to 3.6 X 10^^ cm~^ from the Cvi(i_2) hue and 1.7 x 10^^ cm~^ from Cvi(i_5), which under 
assumption of all carbon being in the form of Cvi ions and with solar abundances translates 
into the lower limit for the warm absorber column equal 1.1 x 10^° cm~^ and 5.1 x 10^^ cm~^, 
correspondingly. 

Exact determination of the column density is complicated. At the flat curve of growth 
we need additional information about velocity dispersion b to indicate column densities. 
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In the case of considered data this is very difficult, because la error of full width of half 

maximum (FWHM) is too high to estimate b parameter. For instance for Sxil line, presented 
in Fig. 3, FWHM = 5ti^ x 10^ km/s. This fact also precludes us to integrate over fitted 
line profile in the aim of determining exact fiux emitted in this line. So, we are unable to 
compute exact optical depth of the line and therefore the column density of absorbing ion. 
The number of lines is to small to fit full model (for instance xstar). 

Another lower limit on column density of Cvi ion can be derived knowing that the 
optical thickness of the Cvi lines is higher than 1. In such a case we can write: 



where T is the temperature of the medium and Ai is atomic mass of the ion. For the Cvi(i_5) 
line we obtain: 



For the temperature approximately equal to 10^ K this value is twice higher than the one 
determined by Mason et al.(2002) in case of Mrk 766. This hmit is over an order of magnitude 
weaker than the limit obtained directly from Equation 1. FWHM estimated from the T ~ 10^ 
K is quite narrow and the line width is rather determined by much larger turbulent motion. 

The value of reduced was equal 0.6363 after fitting all possible absorption fines. 



After absorption lines have been identified, we were able to fit a few emission lines from 
an accretion disk. We applied the standard DISKLINE model with fixed power law index in the 
radial dependence of the irradiating fiux (3 = 2.7. We assumed r^n = 6Rg and Tout = lOOOi?^ 
(where Rg — GM/c^). We have found the best value of inclination angle i = 30° for the fine 
Cv, and we adopted it (keeping frozen) for the other fines. Results are fisted in Table 3, and 
shown in Figs. 1, and 10. 

The last column of Table 3 presents theoretical values of EWs of the same emission lines 
computed for accretion disk illuminated by hard X-ray flux equaled to the local disk flux (full 
model presented in Rozanska et al. 2002 a,b). The model takes into account harder incident 
radiation (F = 1.9), than it is seen in case of Ton S180, and solar metal abundances. For 
moderate illumination in the model we have found for instance the equivalent width of Nex 
equal 11 eV, Ovill equal 18.2, Mgxil equal 10.5, and Fexvil about 8 eV, or many others at 
the same order of magnitude in respect to the reflected continuum (factor of 2-3 smaller in 





Ncvi > 6.33 X 10 



,13 Vf. 



(5) 



3.2.3. Disk emission lines 
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respect to the total continuum). Models do not predict lines with EWs higher than ~ 20 eVs, 
when the disk is in the hydrostatic equilibrium (Rozanska et al. 2002 also figures 

in Nayakshin, Kazanas & Kallman 2000, Ballantyne, Ross & Fabian 2001). Higher values 
can be only obtained from constant density medium, as discussed by Ballantyne, Ross & 
Fabian (2002), particularly when adopting enhanced metalicity. Emission lines also become 
by a factor of a few stronger when the dynamical flare model is adopted, i.e. the time-scale 
of the flare is shorter than the time-scale for a disk surface layer to achieve a hydrostatic 
equilibrium (CoUin et al. 2002). However, in that case lines should be narrow and not 
represented well by a diskline shape. 

Soft disk line detection is therefore a promising tool to constrain the conditions within 
the cold reprocessing medium. Detection of standard broad disk lines presented in this paper, 
if real, suggest no strong departure from hydrostatic equilibrium and the extension of the 
cold disk roughly down to the marginally stable orbit. 

After all lines have been fltted, we flt again broken power-law with Galactic absorption. 
Results presented in Table 1 last verse indicate that fines do not affect continuum very 
strongly. The final value of reduced appeared to be 0.5531. 
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1 

E [keV] (observed) 



Fig. 1. — The LETG spectrum of Ton S180. Solid line denotes total model after broken 
power law, Galactic absorption and all lines have been fitted. 



Table 1. LETG continuum fit in the 0.1-6. keV energy range 



Model NHiGai.) 10'° cm-2 T, E^reak keV r2 xVd-O.f. 



Plain Cent. 1.82771^:^1^ 2.755l|]:[J|7 ^-^^^ o wl ^-^Q^ o oll 1185/1681 

Fixed NniGai.) 1-52 2.689tn?9 l-77llg;}i^ 2.2611^70 1197/1682 

With lines 1.527l^:}f^ 2.6821^:^19 1.6041^:147 2.2951^564 932/1624 
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Fig. 2. — Observed spectrum with the best fit model for 0.19-0.25 keV 
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Fig. 3.— As Fig. 2, but for 0.282-0.32 keV 
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Fig. 4.— As Fig. 2, but for 0.32-0.4 keV 
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Fig. 5.— As Fig. 2, but for 0.4-0.6keV 
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Fig. 6.— As Fig. 2, but for 1.12-1.37 keV 
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Fig. 7.— As Fig. 2, but for 1.5-2 keV 
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Ton S180 3c273 




Energy (keV) Energy (keV) 

Fig. 8. — The comparison of Chandra data of Ton S180 to the Chandra data of 3c273 for 
energies between 0.20 to 0.24 keVs. Similar structure in residuals around 0.22 keV is present 
in both objects. 




Fig. 9. — The comparison of Chandra data of Ton S180 to the Chandra data of 3c273 for 
energies between 0.4 to 0.6 keVs. Two sigma structure presented in Ton S180 around 0.44 
keV is clearly absent in 3c273 data. 
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0.4 0.6 0.8 

E [keV] (observed) 



Fig. 10. — As Fig. 2, but for 0.22-1.05 keV. Three emission lines fitted to the spectrum of 
Ton S180 are seen. 



Table 2. Absorption lines fitted to the Chandra data of Ton S180 



Observed energy (keV) 


EW (eV) 




Best identification 


0.2188in?4 


1.39 ±0.87 


7.42 


instumental 


n 9OOO+0 0011 

U.z>i>iZ_0 0010 


2.05 ± 1.38 


6.6 


instumental 


n 997'O+U-UuUo 

U.ZZ^Z_Q QQQ9 


3.62 ± 1.29 


22.29 


instumental 


0.24141°:^ 


1.87 ± 1.08 


6.34 


instumental 


0.28841^^°^ 


1.82 ± 0.83 


15.35 


instumental 


n OQ9Q+0.0012 
u./:yz,o_o 0012 


1.39 ± 1.17 


6.92 


instumental 


u.zy / '±_o 0022 


6.31 ± 3.70 


14. 


Mgxi (0.2971; 0.3155) 


n "^nvi +0-0026 

U.OU / l_o.o030 


2.26 ± 1.95 


5.18 


Sixiii (0.30878; 0.3279) 


n 47+0-018 

U.O14(_0 030 


2.05 ± 1.28 


7.55 


Sxil (0.31899; 0.3387) 


n Qosc;+o-ooo9 

U.O-iOU_o.oo06 


0.39 ± 0.50 


3.47 


Cv (0.33386; 0.3545) 


U.04Zi_o 0045 


1.38 ± 1.30 


4.45 


Cvi (0.3459; 0.3673) 
blend with Mgxil (0.368) 
Sixiv (0.3703), Cv (0.370£ 


U.44ZU_o.024 


2.22 ± 1.58 


5.16 


Cvi (0.4428; 0.4702) 


0.6963in^^ 


0.74 ± 0.69 


3.47 


Fexvil (0.6868; 0.72932) 
blend with Sxvi (0.731) 


-1 -149+0.0034 
1.14Z_0.oo95 


0.992 ± 1.18 


3.29 


Nex (1.1387; 1.209) 



blend with Fexxv (1.206), 
Fexxv (1.208), and Fexxv (1.213) 
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4. Proposed Model of the Warm Absorber 

In this section, we present models of warm absorber, allowing density, temperature and 
so consequently ionization state to be stratified (Krolik 2002). The only assumption which 
we made is the constant pressure. Such an approach is more realistic than the constant 
density model, since thermal instabilities can produce gas with steep ionization structure 
keeping all zones in pressure equilibrium. 

We have used code titan (Dumont, Abrassart & Collin 2000) to compute temperature 
and density structure of the warm absorber together with radiative transfer. The code 
is based on two-stream approximation and works in plane-parallel geometry. Radiative 
transfer is computed both in lines and in continuum so the code can be used for very 
inhomogeneous even optically thick media. Thermal, ionization and statistical equilibrium 
of ions are computed in complete non-LTE for 10 elements and all corresponding ions. We 
consider following elemental abundances: H: 1, He: 0.085, C: 3.3 x 10^"^, N: 9.1 x 10"^, O: 
6.6 X 10-^ Ne: 8.3 x 10-^ Mg: 2.6 x IQ-^, Si: 3.3 x 10~^ S: 1.6 x IQ-^, Fe: 3.2 x 10"^. We 
transfer totally 721 lines of different ions, taking their oscillator strengths and energies from 
National Institute of Standards and Technology (NIST, http://physics.nist.gov). Photo- 
electric absorption is fully taken into account. 

In our models, we have assumed the shape of incident radiation consistent with this 
fitted to the data of Ton S180. We have adopted power law photon index T = 2.6 (Sec 3.1), 
and spectral limits from 0.012 up to 100 keV as we have concluded from broad band spectral 
energy distribution presented by Turner et al. (2002). 

The assumption of constant pressure refers to the total: radiation and gas pressure, 
and its value is self consistently determined by the code. The ionization parameter on the 
top of the absorber and the total column density are free parameters of our models. We 
have computed models for = 10^ and A^t = 1.7 x lO^^cm"^ named Al, = 1000 and 
Ntot = 6.8 X lO^^cm-^ named A2, and = 10^ and A^t = 1.7 x lO^^cm-^ named A3. Such 
values of column densities for given are upper limits above which optical depth of cloud 
is large enough that strong absorption edges become present in outgoing spectra. 

Even for high , the maximum Compton temperature is low, Tcompt — 2.1 x 10^ K, 
because the incident spectrum is soft (F = 2.6). 

In Fig. 11 we show the temperature structure for Al (dashed line) (higher does not 
change this structure) compared to A2 (solid line) and in Fig. 12 we show the corresponding 
spectra. In case of higher (model Al) the temperature range within the medium is 
much broader and the effective spectrum from such region (presented in Fig. 12 dashed line) 
displays shghtly stronger fines than for A2 (solid line spectrum) . We have enlarged range of 
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Table 2 — Continued 



Observed energy (keV) 


EW (eV) 


Ax' 


Best identification 




1 QOQ+0.0030 
1.00y_o 0024 


1.48 ± 1.84 


2.95 


Mgxil (1.3861; 1.472) 


1 709+0-030 
1.1 Uzi_o 032 


1.88 ± 2.38 


3.7 


Mgxil (1.7327; 1.840) 
blend with Sixill (1.853) 


1.8601°:°!^ 


10.6 ±5.24 


6.09 


Sixiv (1.8831; 2.000) 



*^*^ Expected observed energy and the rest frame energy in keVs for 
identifined hnes are given in parenthesis, for blends only rest frame energy 
is hsted. 



Table 3. Emission disk lines fitted to the Chandra data of Ton S180 



Observed energy (keV) EW (eV) Ax^ Best identification EWs from illuminated disk* 

0.2704l°:°°g7 15.3 ±2.63 72.73 Cv (0.2899; 0.3079) 8.69 

0.6300l^:°J^ 13.5 ± 4.89 21.26 Oviii (0.6149; 0.653) 19.0 

0.94031^21 30.7 ±7.41 46.26 Nex (0.961; 1.020) 11.29 



* Theoretical EWs are computed in respect to refiected continuum 
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frequencies to illustrate region observed in Ton SI 80 data. 

We have also performed comparison of Al model to the warm absorber with constant 
density (Fig. 11 dotted-dashed line) assuming that both clouds have the same total column 
densities and they are illuminated by the equivalent fluxes. Thus, choosing the volume 
density for constant density cloud, the ionization parameter is determined by the value of 
flux the same as for stratified cloud. We have chosen volume density equal 3 x 
which implies .^o = 300. Spectrum for constant density cloud presented in Fig. 12 by dotted- 
dashed line displays strong edges and cannot be applicable to explain our observations. 
Decreasing volume density we increase ionization parameter, and it may remove edges. We 
plan to make proper calculations in the future work to show how to distinguish between low 
density constant cloud and stratified cloud. 

We have also calculated theoretical values of the equivalent widths of transfered lines, 
the same which are found in Ton S180 data. Results are presented in Table 4. EWs for high 
surface ionization (Al) are almost a factor of two higher that for A2 model. In our data we 
see slightly larger EWs and we cannot model them by increasing since model A3 shows 
the same EWs as Al. 

Therefore, we have to introduce model labeled A4 in Table 4, which has exactly the 
same parameters as A2 i.e. .^o = 1000 and Ntot = 6.8 x lO^^cm"^ and additionally we 
assume turbulent velocity Vturb/c = 3 x 10~^. Turbulent velocity does not change the total 
temperature structure, it only affects hues, making their EWs higher, as shown in Table 4. 
For model A4 several lines, for instance Mgxi, Sixill, Sxil Cvi, and Sixiv) have EWs in 
quite good agreement with those achieved from observations (see Table 2). 

In all presented models the ratio of £'l^[Cvi(i_2)] / £'VF[Cvi(i_5)] equals 0.7, in very 
good agreement with observed value of 0.6. 

5. Broad band continuum 

In order to draw attention to the possible constraints resulting from detection of the 
broad disk lines we model the broad band spectrum of Ton S180 assuming three basic spectral 
components: (a) an illuminated accretion disk (b) an optically thick Comptonizing plasma (c) 
an optically thin (non-thermal) Comptonizing plasma. For low or moderate accretion rate we 
take the simplest Keplerian disk model, emitting locally as a black body and parameterized 
fully by the black hole mass, M, and accretion rate, m, given in Eddington units (taking 
efficiency 1/12 appropriate for Newtonian model). We consider also highly supercritical 
accretion disk applying the model described by Kawaguchi (2002). This model includes the 
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Table 4. Equivalent width [eVs] of absorption lines for several 
models computed by titan 



fon 




Transition 


Al 


A2 


A3 


A4 


Sviii 


o 2 o 5 

2s 2p^ 


2"r)0 o 2 o 5 oJ 2t 

P - 2s 2p° 3d L 


0.522 


0.293 


0.506 


1.431 


Fexv 


o 2 Ici 

Ss"^ ^S ■ 


- 3s4p ^P^ 


0.406 


0.191 


0.392 


1 -t A 

1.164 


Six 


28^ 2p 


2ti0 ri 2 o 1 2t^ 

P ^ - 2s 3d D 


0.584 


0.286 


0.568 


1.193 


Sx 


o 2 o 3 

2s^ 2p'' 


4on o 2 o 2/3"r>\ 01 4"r> 

^S" -2s^ 2p^(''P) 31 ^P 


0.585 


0.305 


0.567 


1.549 


Cv 


1 2 1 ri 

Is"' ^S ■ 


- 2p ^P" 


1.112 


0.587 


1.079 


1.849 


Sxi 


2s 2p^ 


3t-> 2 oj 3t 

•^P - 2s 2p 3d 


0.680 


0.333 


0.659 


1.500 


Mgxi 


2p ''P^ 


A 1 3t 

- 41 ''L 


0.834 


0.468 


0.809 


1.995 


Sixiii 


2s - 


o InO 

3p ^P^ 


0.783 


0.439 


0.759 


2.019 


Sxii 


2s"' 2p 


^pu _ 2s^ 3d "'D 


0.719 


0.348 


0.697 


1.580 


Cv 


1S2 IS ■ 


- 3p ipo 


1.287 


0.708 


1.249 


2.151 


Cvi 


Is ^S- 


21 2L 


1.186 


0.599 


1.136 


1.772 


Cvi 


Is ^S- 


51 2l 


1.650 


0.881 


1.585 


2.585 


Fexvii 


2s2 2p6 


^S - 2s2 2p53s ip° 


1.216 


0.621 


1.159 


4.008 


Nex 


Is ^S- 


31 


3.151 


1.856 


3.063 


7.485 


Mgxii 


Is 2S- 


21 2L 


3.561 


2.101 


3.561 


9.107 


Mgxii 


IS 2S- 


41 2L 


5.224 


2.947 


5.330 


12.26 


Sixiv 


Is 2S- 


21 2l 


4.594 


2.426 


4.598 


11.56 



^ 1 means that the configuration is a mixture of multiplets s p d 



-26- 



effect of strong radial advection and the departure from the local black body emission due 
to electron scattering which are essential at very high accretion rates. 

We assume that the energy generated in the inner disk illuminates outer disk, with 
the irradiating flux oc r"^"^, with (5 — 0.1. This irradiation effect may come from the 
scattering of the disk radiation back toward the disk by the warm absorber as well as from 
the direct interception of the radiation by the flaring outer disk. We adopted a simple 
parameterization and did not model this irradiation since the process is quite complex. It 
would require assumption about the radial and angular distribution of the density of the 
scattering medium (see e.g. Kurpiewski et al. 1997) and the computations of the viscosity 
dependent disk vertical structure, as in Rozanska et al. (1999). We did not compute the 
disk vertical structure in order to verify whether the outer disk is flaring strongly enough to 
allow for such a strong effect. 

We assume that the inner disk for (r < rskin) is fully covered by the optically thick 
plasma. Thermal plasma parameters are: the optical depth, r, and the temperature, T. Only 
soft photons from the inner disk come through this plasma, and we compute the spectrum 
using the simple code based on Sunyaev & Titarchuk (1980) description (used for example 
by Collin- Souffrin et al. 1996). We adopt the plane parallel geometry and the generation 
of the photons at the bottom of the zone. This component forms a soft X-ray component 
of the spectrum. We next added a contribution from much hotter, or non-thermal, plasma 
component in order to provide the appropriate hard X-ray tail but the parameters of this 
component are poorly constraint and we do not address this issue in detail. 

Our model provides a one-parameter family of solutions for the broad band spectrum, i.e. 
we can find a good model for any mass of the black hole between 1 x lO^Mo and 1 x lO^Mo, 
but for a fixed mass all other parameters (m, r^fej^, r, and T) are uniquely determined. The 
appropriate accretion rate is needed to reproduce the normalization of the spectrum in the 
optical band, where it is mostly given by the product of M^m ( for illumination case). The 
value of Tskin must be chosen in order not to overproduce/underproduce far-UV flux. The 
slope of the soft X-ray spectrum gives one constraint on the combination of r and T, but 
for low values of T which are always derived the normahzation of the soft Comptonized 
spectrum also depends significantly on the plasma parameters (for a number of soft photons 
already specified by rskin)- 

Examples of the solution are therefore as follows: for M — 1O^M0: we have m — 0.07, 
fskin — '^^Rg, T — 17,T — 0.65 kcV, and for high accretion rate case: M — IO^Mq: we have 
m = 83, Tskin = 140i?y, r = 2.5, T = 12 keV; The broad band data fit for the largest mass 
is shown in Fig. 13 and for the smallest mass — in Fig. 14. 
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Those two extreme models differ considerably with the extension of the disk part covered 
with the hot skin. In the case of small mass, the disk model would overpredict the far-UV 
data points considerably if the large portion of the inner disk is not fully covered with the 
thermal plasma. Plasma parameters are such that they make formation of any emission lines 
rather difficult. 

In the case of high mass the inner disk does not really have to be covered, as indicated 
by much smaller value of Vskin- Clumpy corona model would be easily applicable to this 
case, or even the X-ray emission could come from some shocks in outflowing material at the 
rotation axis. The fraction of the disk covered may be always further reduced assuming that 
synchrotron emission of the plasma itself provides additional soft photons. In such a case 
the lines from the irradiated disk are easily expected. 



6. Discussion 

6.1. Warm absorber 

Reanalyzing Chandra data of Narrow Line Seyfert 1 Galaxy - Ton SI 80 we were able to 
identify several narrow absorption lines due to the warm absorber. Such lines were neither 
discussed by Turner et al. (2001) and Turner et al. (2002) where they presented the data 
for the first time nor seen in XMM-Newton data (Vaughna ct al. 2002) because the lines 
are indeed very faint. We detected lines of highly ionized carbon and other elements, but 
we do not seem to detect any oxygen lines while in the FUSE data (Turner et al. 2002) Ovi 
narrow absorption line components were seen. 

Five of those absorption lines detected in Ton S180 (Cvi(i? = 0.3673 keV), Nex, 
MgXIl(i? = 1.471 keV), MgXU{E = 1.840 keV), and Sixiv) are also seen in another Seyfert 
1 galaxy - NGC 3783 analyzed by Kaspi et al. (2002), with similar values of EWs (0.46, 
2.3, 4.3, 1.6 and 6.6 eV respectively). Also Kaastra et al. (2002) in NGC 5548 have found 
the same Cv, Cvi{E = 0.3673 keV), Nex, Mgxil(£; = 1.471 keV), and Sixiv, with following 
EWs - 0.26, 0.99, 0.47, 3.09, and 3.88 eV. 

We found a lower limit for the warm absorber column of 5.1 x 10^^ cm~^ from the 
equivalent with of the one of C VI lines (see Sect. 3.2.2) and an upper limit of 1 x 1.7 x 10^^ 
cm~^ from our model of the warm absorber (see Sect. 4). The upper limit is not firm since 
it is sensitive to the adopted turbulent velocity as well as abundances. 
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6.2. Broad-band continuum 

The Chandra data show that the effect of the warm absorber does not modify signifi- 
cantly the X-ray spectrum in the case of Ton S180. Therefore our analysis do not change the 
overall picture of the Ton SI 80 continuum, as shown by Turner et al. (2002) and discussed 
by Vaughan et al. (2002). In particular, the huge soft X-ray excess is reasonably explained 
only within the frame of the Comptonization model. 

Variability arguments presented by Vaughan et al. (2002) in favor of hybrid plasma 
model (thermal plasma being responsible for the soft X-ray slope and non-thermal plasma 
component leading to hard X-ray power law) are convincing. However, the data do not 
allow to determine all the hybrid plasma parameters uniquely, as explained by Vaughan 
et al. (2002), so in the present paper we have used much simpler model in order to draw 
attention to the possible constraints resulting from detection of the broad disk lines. 

Those emission lines can be present in the broad band spectrum rather for higher black 
hole mass M — 10^ M© and accretion rate m — 0.07, since with those assumptions hot 
skin does not cover disk completely. This conclusion is opposite to the black hole mass 
estimations done for Ton S180 by Wang & Lu (2001), where M = 1.3 x lO^M©. But our 
result is consistent with studies done by Siemiginowska (1997), where simple disk corona 
model was fitted to the Beppo ^^Xdata (those presented by Comastri et al. 1998). 

7. Conclusions 

We argue, contrary to the previous papers by Turner et al. (2001) and Vaughan et al. 
(2002) that spectral features in absorption and even in emission are present in the Chandra 
data of Ton S180. The quality of the data does not allow us to determine column densities 
of individual metals, we can only put lower limit on column density of Cvi. 

The lack of edges in the spectrum gives upper limit on total column density, which is 
1.7 X 10^^ cm"^ for the best computed model (^o — 10^). 

A few emission lines fitted to the Ton SI 80 Chandra data have lower equivalent widths 
by factor of 6 than those reported by Mason et al. (2002) for Seyfert I galaxy Mrk 766, but 
of similar order as disk fines found by Kaastra et al. (2002) in NGC 5548. 

Estimation of the mass of the black hole in Ton SI 80 based on FWHM of the iJ^j fine and 
the formula of Kaspi et al. (2000) for the relation between the size of the Broad Line Region 
and the object luminosity at 5100 A rather suggests small black hole mass, M = 1.3 x lO^M© 
(Wang & Lu 2001), in contradiction with our conclusion based on the detection of the soft 
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X-ray lines. 
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Fig. 11. — The temperature distribution for models with = 10^ and A^jot — 1-7 X lO^^cm"^ 
(dashed line), ^ = 10^ and Ntot = 6.8 x 10^^ cm^^ (solid line). Constant density case for 
^ = 300 and Ntot — 1.7 x lO^^cm"^ is presented by dotted-dashed line. 
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Fig. 12. — Spectral energy distribution (spectral resolution 60) for models with ^ = 10^ and 
Ntot = 1.7x lO^^cm-2 (dashed line), ^ = 10^ and Ntot = 6.8x lO^^cm-^ (sohd Une). Constant 
density case for ^ = 300 and Nt^t — 1-7 x lO^^cm"^ is presented by dotted-dashed line. Both 
solid and dotted-dashed lines are arbitrarily moved to see differences. Dotted line indicates 
the slope of incident continuum. 
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Fig. 13. — The model of the broad band spectrum of Ton S180 (sohd hne) for the minimum 
acceptable accretion rate. Dotted line gives the spectrum of the outer uncovered part of 
the disk and dashed line shows the radiation from the inner disk after passing through the 
Comptonizing medium. Model parameters: M — IO^Mq, m = 0.07, Tskin — "^ORg, t — 17, 
T = 0.65 keV. 
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Fig. 14. — The model of the broad band spectrum of Ton S180 (sohd hne) for very high 
accretion rate (Kawaguchi 2002) . Dotted hne gives the spectrum of the outer uncovered part 

of the disk and dashed hne shows the radiation from the inner disk after passing through the 
Comptonizing medium. Model parameters: M — IO^Mq, m = 83, Vskin — 140i?p, r = 2.5, 
T = 12 keV. 



